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IY UMa: accretion disc evolution after superoutburst 
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Abstract. CCD photometry of the newly discovered eclipsing dwarf nova IY UMa is used to study the physical 
properties of the accretion disc in the late decline and quiescence stages. Eclipse mapping analysis shows that 
in these stages the accretion disc is cool with approximately flat radial brightness temperature distribution and 
Tbr ~ 5000 — 5500 K. The hot spot is found to lie close along the stream trajectory at distance 0.36a and 0.25a 
from the disc center in February and March. The orbital hump strength decreases by ~ 30% in a month. Taking 
the smallest possible size of the hot spot, we derive an upper limit for its brightness temperature, ~ 15900 K and 
~ 13800 K, respectively. 
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1. Introduction 

An intensive photometric monitoring of the recently dis- 



covered dwarf nova Tmz V85 (Takamizawa 1998) was 
started after January 13, 2000 when Schmeer (200C) re- 
ported a new outburst. The star has been classified as a 
deep eclipsing SU UMa type dwarf nova by Uemura et al. 
( |2000D and got the designation IY UMa (Samus [2000] ) . 
Using numerous long photometric series covering the 
whole superoutburst, Patterson et al. ( 2000| ) determined 
the system parameters: orbital period P Q rb =0*?0739092, 
masses of the two components Mi = 0.68Mq and A/2 = 
O.IOMq, inclination i=85.9° and distance D=190 pc. 
During the superoutburst the star showed both "normal" 
and "late" superhumps with a period of ~0'?07583. 

Here we report our time-resolved photometric observa- 
tions of IY UMa started about the end of the late super- 
humps stage and continued after the accretion disc shrank 
back in quiescence. 



2. Observations and data reduction 

The observations of IY UMa were carried out at Rozhen 
Observatory with the 2.0-m telescope and its attached 
Photometries 1024x1024 CCD camera. The star was ob- 
served in the Johnson V band for three nights, with an 
exposure time of 60 s on February 6 and March 11, 2000 
and with 90 s on March 30, 2000. In addition, B ex- 



posures were regularly taken during the nights. In each 
case, there was a ~15 s dead-time. In total, seven eclipses 
were covered: three on February 6 and March 11 and one 
on March 30. After de-biasing and flat-fielding the pho- 
tometry was done with the standard DAOPHOT proce- 
dures (Stetson 1987). On February 6 and March 30 the 
seeing was ~2 and we performed PSF fitting photome- 
try, which allowed to estimate the magnitudes of the faint 
NE companion: 7=19.5 and B - U=0.69. The data on 
March 11 were processed using aperture photometry in- 
cluding this star and subtracting its flux afterwards. For 
all the observations the stars ID 9 with U=17.633 ±0.037 
and £-7=1.145 ±0.027 and ID 10 with U=17.844 ±0.016 
and B - U=0.620 ±0.136 from the list of Henden (120001) 



served as a comparison star and check, respectively. The 
instrumental b — v colors were transformed to standard 
ones in the usual way with a second-order extinction term 
k b _ v =— 0.035 included. All the B frames were taken out- 
side eclipse, therefore the B — V colors are representative 
for these orbital phases only. The B — V color was found 
to be constant with the orbital phase and averages ~ 0.0 
and ~ ±0.22 for February and March observations, respec- 
tively. These values were used to convert the instrumental 
v band measurements to standard ones. The magnitudes 
were converted to fluxes using the absolute calibration of 
the UBV system (Straizys |l~977|). 
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3. Accretion disc eclipse mapping 

In Fig. Q the light curves are shown, folde d wi th the 
ephemeris determined by Patterson et al. ( 2000| ). The 
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Fig. 1. Photometric light curves, folded with the orbital 
period P or b =0^0739092. The low order polynomial fits to 
the out-of-eclipse data are also shown. 
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analysis of these curves was performed using the system 
parameters determined by Patterson et al. ( 20001 ). To re- 
construct the accretion disc V band brightness distribu- 



tion we used the eclipse mapping method (Home 1985) 
and the Maximum Entropy technique as described by 
Skilling&Bryan ( |1984|) w ith a full azimuthal floating de- 
fault image (Home |1985| ). Since our data were obtained 
with exposure times longer than the WD ingress and 
egress (^35 s, Patterson et al. 200C), we were unable to de- 
compose the light curves into its individual components: 
WD, hot spot and accretion disc. Thus, the curves had 
to be analysed without subtracting the WD and hot spot 
flux. Instead, we accounted for the variation of the hot 
spot flux with the orbital cycle by first fitting the out-of- 
eclipse data with a low-order polynomial (shown with full 
lines in Fig. Q), then normalizing the data to this fit and fi- 
nally rescaling the data to the expected hump flux at zero 
phase. The curve obtained on March 30 was excluded from 
the analysis because of its long exposure time, hence its 
very poor orbital phase resolution. The final curves and 
the fits obtained are shown in Fig.|| and it can be seen 
that the hot spot egress phase changes from ~ 0.09 in 
February to ~ 0.075 in March. The hot spot egress phase 
obtained from the curve not included in the reconstruc- 
tion (taken on 30 March) coincides with that determined 
from the March 11 data. It seems also that in February 
the white dwarf and hot spot are eclipsed simultaneously 
or at least the hot spot is eclipsed very close after the 
white dwarf, while in March the ingress of the two com- 
pact sources seem to be resolved, as was observed for ex- 
ample in Z Cha and OY Car (Wood et al. |l986| , |l989|) . In 
Fig. H are also shown the mid-egress phases of the hot spot 



Fig. 2. Upper panel: Hot spot egress phases. The trend 
after JD2451578 indicates a decrease of accretion disc 
radius. Vertical line marks the end of late superhumps 
era. Lower panel: Light curves of the reconstructed disc 
images. Dotted line is the reconstructed February light 
curve plotted over March data to see the difference clearly. 
Vertical dashed lines show the hot spot mid-egress phase. 



obtained from our data, compared with those derived by 
Patterson et al. ( [2000 ). It can be seen that our results are 
in agreement with a continuous decrease of the hot spot 
mid-egress phase when the system evolved to quiescence. 

The calculated brightness distributions are shown in 
Fig. H Two compact sources are found and they are recog- 
nised as the white dwarf and the hot spot. The stream tra- 
jectories, calculated by integrating the equations of mo- 
tion (Flannery 1975] ) by the fourth-order Runge-Kutta 
method, are also shown in Fig.|^. In both maps the hot 
spot is found close along the trajectory but at different 
distance from the white dwarf, indicating that the accre- 
tion disc has shrank. The radii in February and March, 
determined by the most bright pixel in the hot spot, are 
rd ~ 0.36a and r^ ~ 0.25a, respectively. In the disc re- 
construction from the March data, something like a third 
spot is seen. We believe that this is not a real structure in 
the disc but is an artifact from the reconstruction. 

Having the brightness distribution calculated we can 
estimate the brightness temperature Tbr of the accretion 
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Fig. 3. Reconstructed brightness distribution of the ac- 
cretion disc from February and March data. Arrows mark 
the direction of the maximal hot spot visibility. 

disc assuming black body emission. This is done by solving 
for Tbr the following equation: 



Jv(X)B x (T BR )/XdX 
1 Jv(X)/XdX 1 ' 



(1) 



Here / stands for the emitted flux, -Ba(Tbr) is the black 
body spectral distribution, v(X) is the response of the V 
bandpass and 5 is the pixel area corrected for the fore- 
shortening. In Fig. [| we show the radial distributions of 
the brightness temperature along the accretion disc radius 
together with the expected temperatures for steady accre- 
tion discs for different mass transfer rates. The inner part 
of the disc seems to be well described as a steady accretion 
disc model with M ~ 1O -11 M yr _1 . If this is the case it 
would be in conflict with other results for dwarf novae in 
quiescence where the accretion disc brightness tempera- 
ture radial distribution is found to be much flatter than 
what is expected for steady discs. Most probably, however, 
this is redistributed flux from the white dwarf as a result of 
the algorithm. In the outer part of the disc reconstructed 
from the February data, the brightness temperature dis- 
tribution becomes almost fiat with Tbr — 5500 K. In the 
March data reconstruction this fiat part cannot be seen 
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Fig. 4. Brightness temperature of the disc. Models of 
steady state discs for different mass transfer rates are 
shown also. 

due to the smaller disc radius and the redistribution of the 
white dwarf and the hot spot fluxes. While Tbr ~ 5500 K 
is a typical temperature for dwarf novae in quiescence, 
the observations and numerical studies show that quies- 
cent accretion discs in these systems are most probably 
optically thin. In this case the emission of the disc signif- 
icantly deviates from black body law and the above value 
of Tbr is a very crude estimation of the effective temper- 
ature of the accretion disc. 

The prominent humps seen roughly between phases 
—0.4 and +0.1 are believed to be produced by the hot 
spot. The hot spot light is assumed to be emitted by a 
circular area perpendicular to the orbital plane, and the 



hump is produced by the foreshortening (Wood et al. |1986 , 
1989). These authors supposed that the size of the spot 



is equal to the cross section of the accretion stream (the 
formula derived by Cook& Warner ( |1984| ) ) and used Eq.[l] 
to estimate Tbr of the spot assuming two models of hot 
spot emission: 

i) the hot spot radiates isotropically from both its sides 
with most of the light from nearest to the disc side 
being obscured from our sight, 

ii) a two component model in which the hot spot itself 
radiates anisotropically from its outer side only and 
another isotropic component. 

Since we cannot decompose our light curves, we used this 
method to estimate Tbr of the spot corresponding to the 
anisotropic part in the second model. The flux / was de- 
termined as the difference between maximal flux and the 
flux at phases when the light of the hot spot is not seen 
(in March data the fluxes before and after the hump are 
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significantly different so we took the mean value). The es- 
timations showed that the strength of the hump decreases 
roughly by ~ 30% between February and March, giving 
T BR ~ 15900 K and T BR ~ 13800 K with an error of both 
values ~ 500 K. These values are an upper limit for Tbr 
because the accepted width of the hot spot is a lower limit. 
The actual size of the spot is probably much larger as it 
is produced by the hypersonic impact between the stream 
and disc. 



Wood, J., Home, K., Berriman, G., and Wade, R., 1989, ApJ, 
341, 974 

Uemura, M., Kato, T., Novak, R., Jensen, L.T., Takmizawa, 
K., Schmeer, P., Yamaoka, H., Hensen, A., 2000, IAU Circ. 
7349 



4. Conclusions 

The analysis of the data shows that in spite of the con- 
stant s ystem brightness after February (Fig. 1 in Patterson 
et al. ( |2000 ) and our data) the accretion disc has not 



reached its equilibrium quiescent state. Its radius grad- 
ually decreases as the system evolves to quiescence from 
~ 0.36a in February and ~ 0.25a in March. This is 
in agreement with both, observations and numerical simu- 
lations of dwarf nova outbursts (see review paper of Osaki 
( |l996|) and references therein) . In the same time the tem- 
perature of the disc itself seems to remain constant. The 
light curves obtained on 11 and 30 March are very similar 
without any apparent changes, which suggests that the 
equilibrium state was most probably reached in the be- 
ginning of March. The eclipse mapping shows a hot spot, 
laying along the accretion stream. The strength of the or- 
bital hump decreases by ~ 30%. Assuming that the hump 
is produced by the hot spot and that it emits anisotrop- 
ically from its outer side, we estimate an upper limit of 
the hot spot brightness temperature Tbr ~ 15900 K and 
Tbr — 13800 K in February and March, respectively. Such 
a decrease could be produced by an enhanced mass trans- 
fer rate during outburst or/and by a change of the hot 
spot structure as the angle of impact and densities of the 
disc and stream vary with the accretion disc radius. 
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